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We investigate the quantum non-demolition (QND) measurement of an atomic population based 
on a heterodyne detection and show that the induced back-action allows to prepare both spin- 
squeezed and Dicke states. We use a wavevector formalism to describe the stochastic process of the 
measurement and the associated atomic evolution. Analytical formulas of the atomic distribution 
momenta are derived in the weak coupling regime both for short and long time behavior, and 
they are in good agreement with those obtained by a Monte-Carlo simulation. The experimental 
implementation of the proposed heterodyne detection scheme is discussed. The role played in the 
squeezing process by the spontaneous emission is considered. 
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I. INTRODUCTION 

Generation of entangled collective spin states such as 
spin-squeezed states (SSS) 0, [|| and Dicke states 
is of first importance for quantum information Q and 
quantum metrology (7rlTo|. A non- linear evolution of 
the system is required to create collective entanglement: 
common approaches exploit interactions between the par- 
ticles [ill, IHi ! probing with squeezed light [l3[ or cav- 
ity mediated interaction |9| . The strong non-linearity of 
a measurement determined by the state renormalization 
can also lead to such states 

Measurement induced squeezing of atomic ensembles 
has been performed using either a Mach-Zehnder inter- 
ferometer and simultaneous interaction with several op- 
tical frequencies [Pc| . or cavity transmission modulation 
[l9l [2(jj . Nevertheless, the experimental realization of 
highly engineered atomic states is still challenging and 
a deep understanding of the underlying stochastic pro- 
cess of the measurement is needed. The dynamic of the 
atomic state collapse can be described as a succession of 
partial measurements. Several methods have been used 
to model this collapse process, like stochastic Schrodinger 
equation [3, [2l| , master equation p2j , or wavevector for- 
malism [14j |. 

In this article, we propose to use the freq uency mod- 
ulation (FM) spectroscopy technique |23l426j to generate 
both spin-squeezed and Dicke states. A laser beam is 
phase modulated to produce frequency sidebands; one 
sideband is placed close to an atomic transition and ex- 
periences a phase-shift passing through the atomic sam- 
ple. The detection of the beatnote at the modulation 
frequency allows then to estimate the atomic population 
of the probed state. 

A thcorical analysis of the measurement process shows 
that a heterodyne detection performed on an atomic sam- 
ple allows to prepare new quantum states. The study is 



based on the wavevector formalism previously introduced 
for homodyne detection 11411 and shows similarities with 
the method used in [13, HM to reconstruct the photonic 
state from the results of a QND measurement sequence. 
We study the quantum trajectories of the atomic state 
caused by the repeated interaction with single-photons 
sent in the QND apparatus. In the weak coupling limit 
an analytical expression is derived for the evolution of the 
state variance during the measurement, it completely de- 
scribes the atomic wavefunction collapse. The result is 
in good agreement with a Monte-Carlo simulation of the 
measurement process. 



II. THE MEASUREMENT PROCESS 



The principle of the QND measurement is presented 
in Fig. [TJ An atomic system interacts coherently with 
an optical probe. As a result, the orientation compo- 
nent of the atomic state (such as the atom number in a 
well defined hypcrfme level), becomes weakly entangled 
with the probe photons, and the photodetection grad- 
ually provides the observer with information about the 
spin state. Continuous observation conditionally reduces 
the orientation uncertainty with respect to the QND 
measurement outcome thanks to the measurement quan- 
tum back-action. In any individual measurement trajec- 
tory, the initial wavefunction will thus collapse in a sin- 
gle squeezed state whereas the ensemble averaged uncer- 
tainty will not be reduced, because the measurement out- 
comes over many QND trajectories are distributed within 
the variance of the initial state. Our model will allow to 
study the dynamics of this collapse by considering single 
photon detection events. 
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FIG. 1: Basic scheme of the measurement process: the phase- 
shift induced by the atomic sample is estimated using the 
repeated interaction with single-photons. The back-action on 
the atomic state increases with the precision of the phase 
measurement. 



A. Heterodyne QND detection 

We consider the QND measurement of the population 
difference between two atomic levels using an heterodyne 
detection. The measurement is realised by estimating 
the phase-shift between the two frequency components 
of the detection beam. This configuration is equivalent 
to a frequency space interferometer in which each path 
corresponds to a single frequency mode. In our model, 
we consider photons individually sent into this interfer- 
ometer and reconstruct the output beatnotc photon after 
photon by comparing their arrival time with a time ref- 
erence. 



Model of the apparatus 



tical phase modulator B is used as a spectral beamsplit- 
ter to generate spectrally mode-entangled single-photons 
(see App. [2}. A local oscillator (LO) drives B with fre- 
quency S1/27T. We consider the driving signal as a noise- 
free classical field. After passing through the atoms, the 
two-channel-optical-state is sent to a single-photon de- 
tector which is assumed to be ideal. The detection of a 
photon in the mode (0) is associated with the operator 
Mo, whereas the detection of a photon in the mode (fi) 
is associated with uq. If the coherence length of the op- 
tical source is much larger than the LO wavelength, the 
two modes (0) and (ft) are in a coherent superposition, 
and the spatial overlap of these two modes generates a 
beating at frequency S1/27T. A pulse generator, in phase 
with the LO, delivers regularly spaced pulses (at time 
interval r = 2tt/CI) used as a reference to measure the 
arrival time of the photons. The histogram representa- 
tion of the photon arrival time allows to reconstruct the 
beating signal. 



2. Atomic state 

The atomic ensemble is a collection of a fixed number 
A^ at of two- levels atoms { |a^} , \bi)}, described as fictitious 

1/2-spins {jx,jy,jz } where: 



j1° = ^(Mtel + k-XM), 
3 { v l} = JdoiXftil-lfciXo,!), 
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FIG. 2: Model for the QND based heterodyne detection. 
Single-photons are sent through a spectral beamsplitter B and 
interact non-destructively with the atomic ensemble (C/qnd). 
The detection time of each photon is measured by comparison 
with a time reference signal. 



The model for the QND apparatus is shown in Fig. [5] 
The input optical state is composed of two spectral 
modes. The probe mode (0) at a frequency u>o/2ir is close 
to the frequency of the atomic transition, and will there- 
fore be strongly phase shifted when traveling through the 
atoms. The reference mode (f2) is at a far-off resonance 
frequency (w + fi) /2ir and will be unaffected. A single- 
photon state |1) in the mode (0) is associated to the pho- 
ton annihilation operator an, and the vacuum state |0) in 
the mode (f2) is associated to ac>. A single sideband op- 



The atomic ensemble is characterized by the collective 
spin operators Jk = Y^iJk > k = x,y, z. A collective 
state \n) is an eigenstate of both J z and J 2 = J 2 + J 2 + J 2 , 
and is called a Dicke state. It verifies: 



J 2 In) 



A^at {N, 



2 V 2 

J z \n) = n \n) 



+ i w 



(4) 
(5) 



where — iV at /2 < n < AT at /2. After N p photons detected 
the atomic state is \ip at (N p )) = J2 n c " C^p) \ n )- The ini- 
tial atomic state is a coherent spin state (CSS) polarized 
along J x , which means an average population difference 
(J z ) — 0, and its explicit expression is (29j : 



Cn (0) 



(6) 



From the Moivre-Laplace theorem, if iV at 3> I the coeffi- 
cients c n (0) are well approximated by a gaussian distri- 
bution: 



c n (0) oc exp 



(7) 
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B. Analytical description of the measurement 
process 

1. Scattering matrix of the interferometer 

The spectral beamsplitter B creates a coherent super- 
position of the two input spectral modes (0) and (O), see 
App. El Its action is described by the following matrix: 



B oc 



T 

Re iQt 



R 

Te iQt 



(8) 



where, in analogy with a spatial beamsplitter, R is the 
probability for a photon in the probe mode (0) to be 
in the reference mode (fl) at the output of B, and T is 
the probability for a photon to stay in the carrier mode. 
We neglect any absorption and dispersion in B so that 
R + T = l. 

When luq + Q is far from any atomic transition, the 
reference beam experiences no phase-shift and the QND 
interaction can be written as 



Uqnd (Jz) 







(9) 



<f> is the optical phase-shift induced by a number diference 
of two (n — 1), and depends on the coupling strength 
of the transition at ujq and on the optical density of 
the atomic cloud. Spontaneous emission induced by the 
probe beam is neglected in this analysis; the reduction it 
causes on the coherence of the atomic state is considered 
in Sec. [GO 

The scattering matrix of this system is S (J z ) = 
t^QND (Jz) B and it furnishes the output modes of 
the interferometer once the input ones are given: 
(u (J z ) ,u n (J z )) = S (J z ) (a ,a n )- 



2. Measurement back-action 

A click at the detector corresponds to the annihilition 
of a photon either in mode (0) or (51). After the mea- 
surement the state of the system satisfies the relation: 

|0 ,0n)®|V>at (N p + 1)) 
cx [«o(J») + un(J»)]|lo,O n >®|^.t(JVi,)>. (10) 

The stochastic recurrence relation which defines the 
atomic state evolution after one measurement follows 
from the expression of S ( J z ): 



IV'at (N p + I)) CX 



+i 



VT - Vr] sin 



4>j z nt 



l^at (n p )) , (ii) 



This expression contains a beatnote signal at frequency 
0/27T with a contrast C = 2\jRT. The reconstruction of 
the beatnote photon after photon is achieved by compar- 
ing the arrival time tk of the fc-th photon with the last 
time reference pulse pr, where p is the number of pulses 
counted between tj= and t = tk- We estimate the 
phase-shift tfk = QStk from the time delay 5tk = tk — pr 
between the last pulse of the time reference and the pho- 
ton detection. 

The atomic state after the detection of N p photons is 
obtained applying N p times the relation of Eq. (TTTI) on 
the initial atomic state j^at (N p = 0)). Expressed as a 
superposition of collective states, it results \ip at (Np)) — 
En ^n,, (n) c„ (0) \n), where 

\F Np (n)| 2 oc Y[ [1 +C cos ((f)n-ipk)} (12) 
fe=i 

is the back-action function [40| . 

Eq. (fT2"]) leads straightforwardly to the probability 
P (ip) to measure a phase-shift tp when the (N p + l)-th 
photon is detected: 

p (v) = E i c « ( n p )i 2 ( J + c cos ^ n - v)) ■ ( 13 ) 

n 

Using this expression and the recurrence relation of 
Eq. (fTTj) . we simulated the quantum trajectories followed 
by the atomic state (see App. [C]). Typical results are 
presented in Fig. [3j 

In the weak coupling limit (4>N ai -C 1), it is possi- 

l 1 2 

ble to derive an analytical expression for the Tn p (n) 
function (see App. |B|) : 



| T Np (n) | oc exp -2M 2 N p ( n 2 + 28(pn/<j> 
where 



5 v= l m „w E 5 <p. 



N p /Nt 



> N P 



(14) 



(15) 



is the average mean position over the followed trajectory, 
and 



M 



(16) 



is the measurement strength. This expression of the 
back-action function contains the complete description of 
the atomic state evolution in the weak coupling regime 
and allows to quantitatively study the atomic wavefunc- 
tion collapse, that is the squeezing process. 



III. DYNAMICS OF THE WAVEFUNCTION 
COLLAPSE 



where t is the time at which the photon has been de- 
tected. 



The atomic state evolution is studied by comparing 
the initial atomic state distribution c n (0) of Eq. (J7J with 
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FIG. 3: Trajectories of the atomic distribution momenta obtained with a Monte-Carlo simulation, once set jV a t = 200, <f> = 10 
rad, and R = T = 1/2. (a) Mean position of the atomic population (J z ) as a function of the number of measured photons. For 
a long enough measurement the state converges to a Dicke state, (b) Variance evolution during the measurement sequence. 
The solid line is the average over 1000 trajectories (only 20 are plotted) while the dashed line is the result of the analytical 
calculation of the variance at short time, given by Eq. ()18[) . At short time scale, the evolution is deterministic as predicted 
by the analytical study. At long time scale, the variances become stochastic but remain bounded. In addition, the average 
variance over the trajectories is below the short time behavior, (c) The histogram of the mean positions at the end of the 
trajectories is compared to the initial distribution (solid line). The measurement process satisfies the Born probability rule 
Pn = \(n\^ t {N p = 0)>| 2 . 



Fm (n) ■ Depending on their relative width, two regimes 
appear during the atomic state evolution: the short time 
limit - where the width of the atomic distribution A j| is 
large compare to one and the atomic state is in a super- 
position of many atomic levels - and the long time limit 
- where the distribution is very narrow and the state is 
split on a few levels of amplitude Fn p (n). The boundary 
between these two regimes occurs for A J 2 ~ 1, that is 
for N v ~ M- 2 . 



no + 1/2 n 




FIG. 4: Extreme cases for long time evolution. The mean 
position of the atomic distribution is in the middle of two 
eigenvalues of J z , (a) or centered on an eigenvalue of J z (b). 



A. Short time limit 

In the short time limit, defined by the condition N p <C 

M , LFjVp is broad and \c n (N p )\ is a gaussian 

distribution characterized by the following mean position 
and variance: 

(j,) = -c 2 e* 2 w/<t>, (17) 

AJj = £ 2 AW4, (18) 

where k 2 = M 2 N at N p is the signal-to- noise ratio and 
£ 2 = 1/ (l + n 2 ) is the squeezing factor. For N p > 1, 
the squeezing factor drops below unity and the initial 
CSS collapses into a SSS as a consequence of the mea- 
surement process. The remarkable result is that A J 2 is 
independent of the stochastic parameter Sip at first order 
and is thus deterministic, as shown with the numerical 
simulation in Fig. [3] (b). 



B. Long time limit 

In the long time limit (N p 3> M~ 2 ) the atomic distri- 
bution is very narrow (AJ 2 -C 1) and the atomic state 
c n {N p ) ~ Fm p (n) is spread over a few eigenstates. The 
evolution strongly depends on the distance between the 
mean value (J z ) and the closest eigenvalue of J z . Two 
extreme cases can be considered (Fig. |4j: either (J z ) is 
in the middle of two eigenvalues, which turns out to be 
an unstable equilibrium state of the system, or it is an 
eigenvalue of J z , which corresponds to the stable equi- 
librium state (the stability of the attractors in treated in 

In the first case the distribution is centered between 
two eigenstates of J z , namely |no) and \uq + 1). The 
atomic state is well approximated by a superposition of 
these two states \ip at (Np)) ~ (|™o) + \no + 1)) /v2, since 
the amplitude of the other states decreases exponentially. 
In this case (J z ) = uq + 1/2 and the variance is A J 2 = 
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(V'at (A^) | Jl l^at (^p)) - (^ 2 > 2 = 1/4, which results to be 
the upper bound value in the long time limit. This limit 
clearly appears in the numerical calculations, as shown 
in Fig. 1(b). 

The second case represented in Fig. [4] (b) corresponds 
to an atomic distribution centered on an eigenstate |no) 
of J z . The atomic state is then a superposition of three 
eigenstates 



l^at (Ay> 



E 

=-1,0.1 



Fn p (no + k) | n + k) 



(19) 



From the normalization condition, the symmetry of the 
state around no, and Eq. (fT3|). we have: 



£*=- 



1,0,1 



\Fn p (n Q + k)\ 



\F N (n ±1)1 = IJ'at (n ) I e 



= 1 

2M 2 N„ 



(20) 
(21) 



Since (J z ) = n , it follows AJ 2 = 2 exp (-2M 2 iVp) , 
which is the lower bound for the variance evolution as 
plotted in Fig. 02(b). 
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FIG. 5: Two possible measurement schemes, (a) One side- 
band is in the middle of the two probed states whereas the 
carrier and the other sideband are far from the transition 
(A FF i <C A c ). (b) Each sideband is close to an atomic tran- 
sition, the carrier is in the middle of the two probed states. 



IV. EXPERIMENTAL CONSIDERATIONS 

The main experimental challenge in a QND detection 
is to limit the spontaneous emission and so the optical 
power while reaching the shot-noise limit with a sufficient 
bandwidth. In this paragraph, we show how a heterodyne 
detection allows to reach the shot-noise of a weak probe 
beam and how the existence of two sidebands rejects all 
the frequency independent common dephasing sources. 
We finally discuss the effect of the spontaneous emission 
on the squeezing performances. 

In a realistic setup, different configurations for the 
measurement can be implemented depending on the rel- 
ative position of the carrier and the sidebands with re- 
spect to the atomic transitions. Two examples of possible 
schemes are presented in Fig. 03 where the case of alkaline 
atoms is considered. 



A. Shot-noise limited detection 

The intrinsic advantage of a heterodyne detection is 
that it makes possible to reach the shot-noise for a weak 
probe beam by beating it with a much stronger reference 
beam. In FM spectroscopy, the reference beam is the 
strong carrier and the probe beam is the weak sideband. 
If N c is the number of photons in the carrier, N s that 
one in the sideband, and A^ the photon equivalent noise 
due to the detection electronics, the signal-to-noise ratio 
of the optical detection is: 

SNR cx J NcNs -. (22) 
V N c + N s +N e v ' 

If N c > N s , N e , then SNR ~ y/T^: by using a strong 
reference, the detection can thus be limited to the shot- 
noise of the weak probe. 



Nevertheless, the reference must be far from the atomic 
transition to strongly reduce its absorption. A high mod- 
ulation frequency f2 is thus required, which results to be 
the main technical limitation of the method. Nowadays, 
photodiodes designed for high-speed optical communica- 
tions provided with integrated transimpedance amplifiers 
with gains of a few kilo-ohms, set an indicative upper 
limit of 10 GHz for the achievable bandwidth. 



B. Common mode noise rejection 

All the systems where the phase induced by the atoms 
on a probe beam is measured by comparison with the 
phase of a reference beam require a precise control of 
all the sources of relative phase noise between the two 
beams. The geometrical splitting of the two components, 
like in a Mach-Zehnder interferometer, requires to pre- 
cisely stabilize the relative length between the two paths. 
On the other hand, if the splitting is in frequency and the 
two components are spatially overlapped, all the noise 
sources acting on the optical path length (like acoustic 
noise, optical index fluctuations and mechanical displace- 
ments) are rejected, since they are common mode. In the 
specific case of a two sidebands modulator, where only 
one sideband is close to an atomic resonance, the signal 
detected on the photodiode is the difference between the 
beatnotes of each sideband with the carrier. The differen- 
tial detection leaves only the atomic contribution to the 
optical index, an active stabilisation of the interferometer 
is not required. 

In (l8| . an amplitude modulation method in a Mach- 
Zehnder interferometer is proposed which similarly re- 
jects the common mode noise using two probe beams. 
Nevertheless, it requires a demanding good spatial over- 
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lap, a well balanced power and an active phase control 
of the two probes; all these features are intrinsic in the 
method we propose. 

The proposed scheme could be further improved by in- 
jecting the probe in an optical cavity. The signal-to- noise 
ratio for the same destructivity is increased by a factor 
square-root of the finesse 25|. The injection in another 
cavity resonance can be used to perform a differential 
measurement when studying amplitude fluctuations of a 
DC signal as shown in [IjJ [20| • In our scheme, a differ- 
ential measurement that rejects the phase noise due to 
cavity length fluctuations is realized by placing the car- 
rier in the middle of two cavity resonances and injecting 
the two sidebands in two different resonances. 



C. Spontaneous emission 

The performance of a QND detection is reduced by the 
spontaneous emission induced by the probe. The atomic 
decay causes a decoherence of the atomic state since the 
detection of the emitted photon allows to determine in 
which state the atom has been projected. The pro blem 
of quantifying this effect has been treated in [I3.l3l[ l32j . 

We now analyse the performance of scheme (a) in Fig. [5] 
when the atomic decay is taken into account. The de- 
phasing induced by the atomic cloud on the optical probe 
is given by [Hj|: 



A 



AttA 



(N 1 S 1 +N 2 S 2 ) 



(23) 



where Nf is the population of the \F) state, A the optical 
wavelength, A the probe beam area, and Sf the probe 
coupling to the \F) state: 

jAff' 



S f = 



/-^ A 2 4- ■y 2 
„, ^FF' + 7 



Sff' 



(24) 



In the last equation 7 is the atomic linewidth and Sff 1 
is the dipole transition strength of \F) — > \F') given by: 



S FF - = (2F' + 1) (2 J + 1) 



J J' 1 
F' F I 



(25) 



where I and J are the nuclear and total angular mo- 
mentum, respectively. The detunings A.if> and A 2 p> 
are choosen so that Si — — S2 = S. Introducing e 
so that Ni = N at (1 + e) /2 and N 2 = A at (1 - e) /2, 
the optical dephasing can be written as A$ = poSe, 
where p = \ 2 N at /4:TrA is the resonant optical density. 
The phase-shift cj> induced by a population difference of 
one atom between the two atomic states is obtained for 
e = 2/N at : 

*"53- < 26 > 

The integrated probability for an atom to scatter a pho- 
ton when probed by N p photons is, in the limit C< 1: 

Po C 2 N P 



V 



A/a 



(27) 



where C is the atomic lineshape: 

,2 



7 



F,F' 



A 2 4- -Y 2 
Ly FF' ' 7 



Sff> 



(28) 



The signal-to-noise ratio can be rewritten as a function 
of the decoherence: n 2 = pp$r\ where (1 = S 2 /C. In [3lj |. 
the squeezing factor is estimated to be: 



e = ^-+i-d-v) 2 . 

1 + PPoV 



(29) 



Considering 10 7 atoms of 87 Rb optically trapped with 
a beam waist of 20 pm and probed on the D 2 (5 2 Si/ 2 — > 
5 2 P 3 / 2 ) line at A = 780 nm with a detuning A13 = 3.2 
GHz then // ~ 1 and = 4.1 10~ 7 rad. The result- 
ing resonant optical density is po ~ 2400 from which a 
squeezing factor ^ 2 = 0.06 (about 12 dB) can be achieved 
for an optimum decoherence 77 ~ 0.01. It corresponds to 
a measurement of N p ~ 10 7 photons in the total beam 
for a modulation depth of 1%. 

The optimum squeezing factor strongly depends on 
the optical density of the atomic sample. Using a Bose- 
Einstein condensate (BEC), a much higher squeezing can 
be achieved but much less photons have to be used during 
the measurement to reach the optimum. 

With regard to the possibility of exploring the long 
time regime behavior to prepare Dicke states, the de- 
coherence due to spontaneous emission can be seen as 
a constraint that would make this regime hard to reach. 
Nevertheless the discussion above holds mainly for a large 
atom number, if the initial sample contains only a few 
atoms (for example a BEC of 100 atoms) the variance 
of the coherent state is of the order of a few units so 
that entering the Dicke state regime in the presence of 
spontaneous emission is realistic. 



V. CONCLUSION 

In this article we analyse a QND heterodyne detection 
of an atomic sample to prepare spin-squeezed states. A 
theorical model of the measurement process is proposed 
that consists on the sequential detection of single pho- 
tons that causes a back-action on the atomic wavevector. 
One should note that in Eq. (TT21 no asumption has been 
done on the atom-light coupling strength. An analytical 
solution is obtained in the weak coupling regime and com- 
pletely describes the atomic state evolution. Remarkably, 
two regimes occur during the wavefunction collapse: at 
short time spin-squeezed states are prepared whereas it 
converges to Dicke states at long time. 

Finally, experimental considerations to implement the 
detection scheme have been presented, pointing out the 
main advantages and limitations of the method as well 
as the role played by the atomic decay and its influ- 
ence on the squeezing factor. The requirement of an 
active stabilisation of the separated geometrical paths 
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used in standard interferometric methods, such as the 
Mach-Zehnder configuration, is replaced by the need for 
a stable RF field. Moreover, the heterodyne detection 
allows to reach the shot-noise of a weak beam by mixing 
it with a stronger one. 

This detection method could be used to generate other 
kinds of non-classical states such as Schrodinger cats or 
NOON states. Such states could be employed in quantum 
enhanced metrology to approach the Heisenberg limit 
sensitivity 34]- Moreover, sequential measurements al- 
low to prepare almost deterministic states by implement- 
ing a quantum feedback loop 0, 135U37I . and can be ex- 
ploited for quantum error correction 38] . 
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The phase modulator thus splits the input photon wave- 
function in two spectral components, whose superposi- 
tion determines a beating of the photon wavefunction 
with itself at the modulation frequency of the RF field, 
as shown in Fig [6l 




Phase 
modulator 



FIG. 6: Quantum beats of the single-photon wavefunction 
generated by a phase modulator. 



Appendix B: Derivation of an analytical expression 
of j^iVp (™)| 2 in the weak coupling limit 



Appendix A: The phase modulator as a 
single-photon spectral beamsplitter 

The electro-optic modulation relies on the Pockels ef- 
fect, a second order non-linear interaction between an op- 
tical and a quasi-static, radio-frequency (RF) field. The 
Hamiltonian associated with this process is: 



The general stochastic measurement process described 
by Eq. (|12[) is non-Markovian. In the weak coupling limit 
(4>N a t < 1) the process becomes Markovian [39( but is 
not stationary. Using again the weak coupling approxi- 
mation we show that the overall process can be split into 
a sequence of stationary sub-processes, which brings to 

I 1 2 

an analytical expression of Tn v (n) . 



Hp = ihg (a^ajai - a m a a[ 



(Al) 



where a m , ao and a\ are the annihilition operators for 
the RF field, the input field and the generated field, re- 
spectively and g is a coupling constant (depending on 
the non- linear crystal, the field frequencies, . . . ). This 
Hamiltonian is the result of two processes: a down con- 
version from an optical to a quasi-static electric field, 
that is an optical rectification, and a frequency sum with 
a quasi-static electric field, which is the Pockels effect. 

Let t be the propagation time of the optical wave in 
the non-linear crystal. The evolution of an incident state 
|^i n ) through the crystal in the weak coupling limit (gt <C 
1) gives the following output state: 



-iHpt/h 



1 



9t (at n al 



(A2) 
IVU ■ (A3) 



Considering a coherent RF field and a single-photon op- 
tical state in input \ipin) — \<X-m) ® | lo, Oi), the output 
field is: 



®(|lo,0i)+s*a m |0o,li». 



(A4) 



1. Evaluation of the probability P (ip) 

In the weak coupling regime the phase shift determined 
by the whole atomic sample is small, then <jm <C 1 for 
-iVat/2 < n < iV at /2. We can expand Eq. (JT31) to 
the first order in <pn around a centered probability dis- 
tribution P ((p) = (1 +Ccos<p) /2tt. We find P(tp) = 
P {(f) +SP [Sp\ {ip), where SP [5<p\ (<p) = -C sin^ <5^/27r 
and Sep — —<f> (J z ) is a stochastic parameter depending 
on the followed trajectory. The chain of stochastic vari- 
ables {<?fc}i<fc<jv * s reduced to a single parameter 5tpN p 
which proves that the process is Markovian. Pq (ip) is 
the deterministic contribution for a zero phase detection, 
whereas 5P [Sip] (ip) describes the atomic phase-shift. 



2. Decomposition into stationary sub-processes 

To study the evolution of the atomic state, we split 
the measurement into sequences of N t photon detec- 
tions, each inducing a negligible evolution of the atomic 
wavefunction (the sub-processes are stationary). A 
small signal-to-noise ratio is required for each sequence 
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{<j) 2 N a tNt < 1). In the weak coupling limit the number where 
of photons in each sequence can be large since (j) 2 N at <C 1. 

The detection probability at the end of the j-th 
sequence of N t measurements is Pj {tp) = Pq (tp) + 
6P[8tpj] (ip). If we choose a phase measurement reso- 
lution ir/m (m € N and m 3> 1) we find from Eq. fj 12[) : 



N p /Nt m 

|^r,(n)| 2 « 1] II f™( n )> ( B1 ) 

j—1 l=—m 



having defined: 



fw' l) ( n ) = I 1 + C cos (<M - 7r//m)] 



(B2) 



Nj i = —Pj(—)N t is the number of measurements giv- 
ing the result tp — irl/m during the j-th sequence. 

The number of photons has to be large in the intervals 
[— ir, — ir + ir/m\ and [it — ir/m, it], where P (tp) is close 
to zero. This is verified for Nt 3> m 3 , since the following 
condition must be satisfied when C ~ 1: 



N 



P {ip) dtp ~ N t 



l-C TT 2 C 



2m 



12m 3 



> 1. (B3) 



In that case Nj t i 1 and the second order term for the 
Taylor expansion of /Jy (n) in 4>n can be identified with 
that of a gaussian distribution: 



(n) ex exp \-2MfN hl (n - m) 



(B4) 



C(/) 2 C + cos (irl/m) 
4 (1 +Ccos(7rZ/m)) 
I 1 + C cos {irl/m 



2 ' 



C + cos {irl/m) 



sin {irl/m) 



(B5) 
(B6) 



By converting the discrete sum over I into an integral and 
using Eq. (jBTJ and Eq. (|B4]). we obtain Eq. (TT4")) . 



Appendix C: Numerical simulations 

The numerical simulation of the quantum trajectories 
shown in Fig. [3] adopts the state of Eq. ([7} as initial one. 
For each photon, the probability P^ p+ i {tp) to measure a 
phase tp for the next detected photon is given by Eq. f| 13[) . 
More explicitly, we use the cumulative distribution func- 
tion associated to this density of probability, which is: 



(¥>) = 



Pi 



N„ 



(6) d6. 



(CI) 



By generating a random number with an uniform distri- 
bution over [0, 1] and numerically inverting fV +i we get 
the phase detected for the {N p + l)-th photon. The re- 
currence relation given by Eq. yields the new atomic 
distribution. The quantum trajectory is obtained iterat- 
ing the sequence. 
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